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Solid angle subtended by a cylindrical detector 
at a point source in terms of elliptic integrals 
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Abstract 

The solid angle subtended by a right circular cylinder at a point source located 
at an arbitrary position generally consists of a sum of two terms: that defined by 
the cylindrical surface (^cyi) and the other by either of the end circles (^drc)- We 
derive an expression for Qcyi in terms of elliptic integrals of the first and third kinds 
and give similar expressions for flcirc using integrals of the first and second kinds. 
These latter can be used alternatively to an expression also in terms of elliptic 
integrals, due to Philip A. Macklin and included as a footnote in Masket (Rev. Sci. 
Instr., 28 (3), 191-197, 1957). The solid angle subtended by the whole cylinder when 
the source is located at an arbitrary location can then be calculated using elliptic 
integrals. 
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1 Introduction 


The knowledge of the solid angle (hi) snbtended by a right, hnite, circnlar 
cylinder at a point isotropic sonrce is reqnired in nnmerons problems in nn- 
clear and radiation physics. Generally, G can be expressed as snm of two 
components: that snbtended by the cylindrical snrface {flcyi) and the other by 
either of the end circles (Gcirc)- Throngh the years this calcnlation has been 
addressed by varions anthors using different methods. Without the worry of 
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being exhaustive we give some examples of such works. iMasketl (jlQhTll out¬ 
lined a general procedure based on Stokes theorem to reduce the double in¬ 
tegral = ff sin 6d6dip to a contour integral in a single variable {6 or ip). 
The method was used to express flcirc and flcyi as single integrals which were 
numerically integrated. Extensive tables describing these result s both for the 
disc and the entire cylinder were reported in a separate work fjMasket et al . 
195fill . An approximati on to the solid a ngle dehned by two parallel plane sur¬ 
faces was described by Gillespie ( IQTlil and applied in the cases of two equal 
rectangles and two equal circles in a face-to-face geometry. With this method 
each surface is subdivided into small finite areas and the two double integrals 
are then replaced by a double sum mation. The calculation of Qcirc was treated 
by Gardner and Verghesel (1971) by replacing the disc with a regular n-side 
polygon of equal area, for which an analytical expression was given. In a sim¬ 
ilar way, flcyi was approxi mated by the analyt i cal expression for an n-side 
regular polyhedral surface ( Verghese et al . Il972ll . iGreen et aA ()l974ll used the 
Monte Carlo method to calculate the cylinder solid angle for two height-to- 
radius ratios (1:1, 2:1), considering distances from source to cylinder center 
up to 12 cylinder radii and angular positions of the source ranging from 0° to 
90" from the cylinder axis. 


An analytic al expres s ion fo r fldrc in terms of elliptic integral s due to Philip 
A. Macklin (jMacklinl. Il957rl appears included as a footnote in iMasketl ()l957ll . 
In the present work we show that also flcyi can be reduced to elliptic integrals 
and give, without derivation, expressions for fldrc which can be deduced in 
a akin way and are different from that due to P.A. Macklin. The solid an¬ 
gle of the whole cylinder can then be expre ssed in terms of ellipti c integrals 
which are rather well known functions (e.g. iMilne-Thomsonl. 119641 for which 
computation algorithms and tables are readily available. 


Instead of blind numerical integration, one can turn to the variety of numeric 
methods already existent, which enable the fast calculation of the solid an¬ 
gle for the whole range of parameters. For instance, the complete integrals 
of the first and seco nd kinds c an be computed using fi) the polynoni i al ap- 
proximations due to iHastingj (|l955ll and included in iMilne-Thomsc^ ()l964 
eqs. 17.3.33 to 17.3.36 ): (ii) the process of the arit hmetic-geometric mean 
f Milne-Thomsonl . Il964 17.6) or (iii) the infinite series ( Milne-Thomson . 1964 . 
17.3.11, 17.3.12) which can be used in combination with Landen’s transfor¬ 
mation when the modular angle is close to 7t/2. 


Since the solid angle can be decomposed into elliptic integrals, any possibility 
of Ending general analytical expressions in terms of elementary functions is 
henceforward precluded. On the other hand, the calculation has been put 
under the sound roof of the subject of elliptic integrals and functions. 


In a recent work ( Prata . 200211 we describe the analytical calculation of the 
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solid angle subtended by a cylinder at a point cosine source. Combining these 
results with those presented here, the case of an axially symmetric point source 
with an angular distribution given by fk{^) = 1 + a 17 ■ k can be treated 
analitically when the source axis (k) is orthogonal to that of the cylinder. 


2 Solid Angle Calculation 


The solid angle subtended by a given surface at a point isotropic source can 
be dehned as 

^surf JJ dVt , (1) 

directions 
hitting surface 

so that 0 < Qsurf < 1. As previously said, in the case of a circular cylinder, the 
solid angle 17 is in general given by 17 = Qcyi + ^drc, where Qcyi and Qcirc are 
the contributions of the cylindrical surface and of one of the end circles. These 
quantities can be dehned with the help of hgs. 1 and 2, where the source is 
located at the origin of the coordinate system and the axis is both parallel to 
the cylinder axis and orthogonal to the planes of the circles. To calculate Qcyi 
it is sufficient to consider the situation depicted in £g. 1, where the source lies 
in one of the planes delimiting the cylindrical surface. Let Qcyio{L, r, d) denote 
the solid angle in this case and flcirc{L, r, d) the solid angle defined by the disc. 
Let us illustrate how 12 can be calculated in an arbitrary situation. Referring 
again to kg. 1, but assuming that the source is located at (0,0, — |; 2 |), one 
has 12 = Qcyio{L + \z\ ,r,d) — flcyio{\z\ ,r,d) + VLcirc{\z\ ,r,d). If the source is 
located at (0, 0, z) and 0 < z < L, then 12 = Qcyio{z, r, d) + Qcyio{L — z, r, d). 
li r > d and z < 0 or z > L, the solid angle is defined by one of the circles 
alone, e.g. 12 = Qcirc{\z\ , r, d) for z < 0. 



Fig. 1. Notation for ^cyio 
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Fig. 2. Notation for Q.circ 

2.1 Calculation of Qcyio 



Fig. 3. Quantities used to calculate 
From figs. 1 and 3, it follows that 

V^max ^max 

^cyio{L, r, d) = 1/(47r) / J sin 9 d9d(p = 1/(27r) / (cos - cos 6'max) dip 

‘Pmin fimin 0 

where 


cos 


cos ^max 
and 


V^max V^min 

L/^L"^ + pf{ip) , 

0 


<~Po 


arcsin(r/(i) , 


( 2 ) 


P-{ip) = dcos ip — 



{d sin ip) 


2 


(3) 
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Thus, 


V^o 

Vtcyio = L/{2T^) j d(p . (4) 

0 

Then we change the integration variable to 7 _ represented in £g. 3 and given 
by 

7- = 7r/2 - 0_/2 (5) 

where 

0-/2 = arctan [sin(99)p _/{r + d — cos((p)p_)] 
and p-(<p) is obtained from eq. 3. 

Eq. 4 is rewritten as 

^cyio = L! (27r) [ = 2/ ^ 1 ^7- (6) 

I ^L^ + pUj.) [p-M J 

where 



(7) 


Introducing 


7o = (7r/2 + ipo)/2 . 


m = 4rd/{L^ + {d + r)^) , 


( 8 ) 


(9) 


n = 4rd/(d + r)^ , (10) 

there results from eq. 7: 

P-i'l-) = {d + r)^ (1 — nsin^7_) , 
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\Jl‘^ + P-(7-) = + {d + rf y 1 - msin^ 7_ , 

L/^ + {d + r)^ = Y^l — m/n 
and 

(d — r)/(d + r) = Vl — n ;d> r . 

Substituting in the rhs of eq. 6 yields 

^cyio = l/(2vr)^l - m/n / . \ ^ 

7o Vi-’"®™ '>'■ 

The integral is easily decomposed into elliptic integrals in the Legendre form: 


y/l — n 

(1—nsin^ 7_) 




^c!,to = l/(27r)\/l - m/n{Vl - n[7T{n; m) - 7r(Tz;7o|m)] - [A’(m) - A'( 7 o|m)]} 

( 11 ) 

where F( 7 o|m), 7 r(n; 7 o|m) are the incomplete elliptic integrals of the hrst and 
third kinds, respectively; and K{m) = F('7r/2|m'), 7r(n-m) = 7r(n;7r/2|m) are 
the corresponding complete integrals (|Milne- Thomson . 1964 1. From eqs. 9 and 
10 it is clear that the parameter m and the characteristic n verify 


0<m<l ,(m = l77L = 0A(i = r) , (12) 


0<n<l ,(n = l77(i = r) 


(13) 


and 


0<m<n<l ;(m = n77L = 0) . 


(14) 


For computational purposes it is worth mentioning that, because of eq. 14, 
both 7r(n:7nlm) and 7r( n: m) in the rhs of eq. 11 fall into the the circular case 
( Milne-Thomsonl . Il964 17.7.11), except when L = 0 or d = r. These special 
cases are treated separately in section 2.3. 


2.1.1 ^cyio when d^r 

From eqs. 8 and 2 results that, as d —> r+, 7 ,, —>• 7 r /2 and the two subtractions 
of elliptic integrals in the rhs of eq. 11 approach zero (e.g. K{m) — F( 7 o|m) —)■ 
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0). T o avoid the consequent loss of precision, an approximation^ (jMasket 
19571 eq. 32) can be used: 


Vtcyio = (27r) ^{(po - l/2[rdcos(po - r^(7r/2 - ipo)]L ^ 

+ 3/8[rd((i^ + 2r^) cosipo — r^(r^ + 2(f‘){'K/2 — ipo)]L~^ “ ■ • 5 (15) 

where the cancellations can be more easily dealt with. Eq. 15 results from 
expanding and termwise integrating the integrand on the rhs of eq. 4. Since 
d — r < p-{‘p) < V d"^ — when 0 < (p < p>o, the expansion can safely be 
done for \/d? — r'^ < L and eq. 15 gives a good approximation (better than 
four digits) for y/d"^ — < L/10. 


2.2 Expressions for Vtcirc 


In a similar fashion to that used to obtain Vtcyio one can show that the solid 
angle subtended by a circular disc (£g. 2) can be written using complete elliptic 
integrals of the third and first kinds: 


Vtcirc{L,r^d) = {2 ti) — m/n{-\/l — n 7r(?7,;7r/2|m) — K{rn)} ]d> r (16) 


Q.circ{L, r, d) = 1/2 — (27r) — m/n{\/l — n 7r(n; 7r/2|m) + K{m)} ]d < r 

(17) 


where m, n are obtained from eqs. 9, 10. 


The complete integral of the third kind can be expressed in terms of incomplete 
integrals of the hrst and sec ond kinds which are e asier to calculate. Since 
0 < m < n < 1 one can use ( Milne-Thomson . 1964 . 17.7.14,17.4.40) to show 
that: 


Etcirc = 1/4 — (27r) ^n/{l + yjl — n)\Jl — m/nK{m) 

— {2'K)~^{[E{m) — K{m)]F{e\rn!) + K{m)E{e\rn!)} ,d>r, (18) 


^ We include here a normalizing factor of l/(47r). 
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^circ = 1/4 — (27r) ^(1 + a/1 — n)^l — m/nK{m) 

+ {2tt)~^{[E{ m) — K{m)\F{e\m') + K{m)E{e\m')} ,d < r , (19) 


where 


m = 1 — m , 


( 20 ) 


e 


arcsin \J{1 — n) /{1 — m) , 


( 21 ) 


E{ip\m) is the elliptic integral of the second kind and E{m) = E{TT/2\m) is 
the complete integral. 


As ail alternative to eqs. 18, 19, one can use a n exp ression ^ due to iMacklin 
( 195711 which appears as a footnote in iMasketl (jl957ll : 


{4:TT)Vtcirc{Li A d) = 27r + 2[K(rn) — E(rn)] [E{6\m) + F('0|m')] 
-2K{m) \^E{6\m') + E('0|m') + 2/l/([l + (a + /9)^]5[/3 + (1 + a^)^])} , 

( 22 ) 

where 

9 = arcsin [^+("+d) ^ = arcsin 

/3+(l+a2)^ [l+(a-/3)2]7 

0 < d < 7r/2 , —7r/2 <'ip< 7r/2 , 0 > \'ip\ , 

0<a = d/L<oo, 0<l3 = r/L<oo 
and m, ml are obtained from eqs. 9 and 20. 


2.3 Special values and continuity 


Since K{m = 1) = oo and 7r(n;7r/2|m) = cx) whenever n = 1 or m = 1, the 
cases L = 0 and d = r have to be studied separately (see eqs. 12 and 13). 

^ We introduce a factor of 47r to conform to the normalization used in the present 
work. 
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Starting from eq. 4 it is straightforward that flcyio{d —r+,L ^ 0) = 1/4. 
Since Qcyio{d > r, L = 0) =0, flcyio is di scontinuous when d — ,L —>• 0. 
Using the properties of elliptic integrals ( Milne- Thomsoni 1964ll one can easily 
show from eqs. 18 and 19 that: 


^circid —^ L 7^ 0) = 1/4 - (27r) Vl - mi K{mi) , 


(23) 


Qcirc{L —> 0) = {0 (d > r); 1/4 (d = r); 1/2 (d < r)} 


where mi = m\r=d = 4r‘^/(L‘^ ^circ is thus continuous except for L = 0. 

Eq. 23 is the same as iMasketl ()l957L eq. 24), divided by a normalizing factor 
of dvr. 


3 Summary and outlook 


An expression for fl^yi in terms of complete and incomplete elliptic integrals of 
the hrst and third kinds has been derived in the case where the source lies in the 
same plane as one of the end discs (eq. 11 ). Expressions for fldrc decomposed 
into complete and incomplete elliptic integrals of the hrst and second kinds 
using a single amplitude (e) were given (eqs. 18 and 19). T hese expressions 
can be used alternatively to one (eq.22) due to iMacklin ( 1957 1 containing also 
complete and incomplete elliptic integrals of the hrst and second kinds but 
using two ampli tudes (9,'ll)). Computat ional methods for the elliptic integrals 
can be found in Milne-ThomsonI (jl964 h The required third kind integrals all 
belong to the circular case (0 < m < u < 1) except for a few cases which 
are treated in section 2.3. The solid angle when the source is at an arbitrary 
position can be calculated using the expressions referred, as described in the 
beginning of section 2. 


The solid angle dehned by a right circular cylinder and an isotropic source 
distributed on a wire parallel to the cylinder axis can also be decomposed into 
elliptic integrals in a way akin to that used here to obtain Qcyio- A work where 
we report this result is in preparation. 
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